Single-Cell Transcription Mapping of Murine and Human Mammary Organoids Responses to Female Hormones.

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Single-cell transcription mapping of murine and human mammary organoids responses to female hormones.

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Clinical applications of 3D normal and breast cancer organoids: A review of concepts and methods.

While mouse models and two-dimensional (2D) cell culture systems have dominated as research tools for cancer biology, three-dimensional (3D) cultures have gained traction as a new approach that retains features of in vivo biology within an in vitro system. Over time, 3D culture systems have evolved from spheroids and tumorspheres to organoids, and by doing so, they have become more complex and representative of original tissue. Such technological improvements have mostly benefited the study of heterogeneous solid tumors, like those found in breast cancer (BC), by providing an attractive avenue for scalable drug testing and biobank generation. Experimentally, organoids have been used in the BC field to dissect mechanisms related to cellular invasion and metastasis-and through co-culture methods-epithelial interactions with stromal and immune cells. In addition, organoid studies of wild-type mouse models and healthy donor samples have provided insight into the basic developmental cellular and molecular biology of the mammary gland, which may inform one's understanding of the initial stages of cancer development and progression.

Innate and Adaptive Immunity to Transfused Allogeneic RBCs in Mice Requires MyD88.

RBC transfusion therapy is essential for the treatment of anemia. A serious complication of transfusion is the development of non-ABO alloantibodies to polymorphic RBC Ags; yet, mechanisms of alloantibody formation remain unclear. Storage of mouse RBCs before transfusion increases RBC immunogenicity through an unknown mechanism. We previously reported that sterile, stored mouse RBCs activate splenic dendritic cells (DCs), which are required for alloimmunization. Here we transfused mice with allogeneic RBCs to test whether stored RBCs activate pattern recognition receptors (PRRs) on recipient DCs to induce adaptive immunity. TLRs are a class of PRRs that regulate DC activation, which signal through two adapter molecules: MyD88 and TRIF. We show that the inflammatory cytokine response, DC activation and migration, and the subsequent alloantibody response to transfused RBCs require MyD88 but not TRIF, suggesting that a restricted set of PRRs are responsible for sensing RBCs and triggering alloimmunization.

Ensuring scientific reproducibility in bio-macromolecular modeling via extensive, automated benchmarks.

Each year vast international resources are wasted on irreproducible research. The scientific community has been slow to adopt standard software engineering practices, despite the increases in high-dimensional data, complexities of workflows, and computational environments. Here we show how scientific software applications can be created in a reproducible manner when simple design goals for reproducibility are met. We describe the implementation of a test server framework and 40 scientific benchmarks, covering numerous applications in Rosetta bio-macromolecular modeling. High performance computing cluster integration allows these benchmarks to run continuously and automatically. Detailed protocol captures are useful for developers and users of Rosetta and other macromolecular modeling tools. The framework and design concepts presented here are valuable for developers and users of any type of scientific software and for the scientific community to create reproducible methods. Specific examples highlight the utility of this framework, and the comprehensive documentation illustrates the ease of adding new tests in a matter of hours.

Prediction of Protein Mutational Free Energy: Benchmark and Sampling Improvements Increase Classification Accuracy.

Software to predict the change in protein stability upon point mutation is a valuable tool for a number of biotechnological and scientific problems. To facilitate the development of such software and provide easy access to the available experimental data, the ProTherm database was created. Biases in the methods and types of information collected has led to disparity in the types of mutations for which experimental data is available. For example, mutations to alanine are hugely overrepresented whereas those involving charged residues, especially from one charged residue to another, are underrepresented. ProTherm subsets created as benchmark sets that do not account for this often underrepresent tense certain mutational types. This issue introduces systematic biases into previously published protocols' ability to accurately predict the change in folding energy on these classes of mutations. To resolve this issue, we have generated a new benchmark set with these problems corrected. We have then used the benchmark set to test a number of improvements to the point mutation energetics tools in the Rosetta software suite.

Macromolecular modeling and design in Rosetta: recent methods and frameworks.

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.

Better together: Elements of successful scientific software development in a distributed collaborative community.

Many scientific disciplines rely on computational methods for data analysis, model generation, and prediction. Implementing these methods is often accomplished by researchers with domain expertise but without formal training in software engineering or computer science. This arrangement has led to underappreciation of sustainability and maintainability of scientific software tools developed in academic environments. Some software tools have avoided this fate, including the scientific library Rosetta. We use this software and its community as a case study to show how modern software development can be accomplished successfully, irrespective of subject area. Rosetta is one of the largest software suites for macromolecular modeling, with 3.1 million lines of code and many state-of-the-art applications. Since the mid 1990s, the software has been developed collaboratively by the RosettaCommons, a community of academics from over 60 institutions worldwide with diverse backgrounds including chemistry, biology, physiology, physics, engineering, mathematics, and computer science. Developing this software suite has provided us with more than two decades of experience in how to effectively develop advanced scientific software in a global community with hundreds of contributors. Here we illustrate the functioning of this development community by addressing technical aspects (like version control, testing, and maintenance), community-building strategies, diversity efforts, software dissemination, and user support. We demonstrate how modern computational research can thrive in a distributed collaborative community. The practices described here are independent of subject area and can be readily adopted by other software development communities.

Structure and function of the immune system in the spleen.

The spleen is the largest secondary lymphoid organ in the body and, as such, hosts a wide range of immunologic functions alongside its roles in hematopoiesis and red blood cell clearance. The physical organization of the spleen allows it to filter blood of pathogens and abnormal cells and facilitate low-probability interactions between antigen-presenting cells (APCs) and cognate lymphocytes. APCs specific to the spleen regulate the T and B cell response to these antigenic targets in the blood. This review will focus on cell types, cell organization, and immunologic functions specific to the spleen and how these affect initiation of adaptive immunity to systemic blood-borne antigens. Potential differences in structure and function between mouse and human spleen will also be discussed.

Loss of Protocadherin-12 Leads to Diencephalic-Mesencephalic Junction Dysplasia Syndrome.

OBJECTIVE: To identify causes of the autosomal-recessive malformation, diencephalic-mesencephalic junction dysplasia (DMJD) syndrome. METHODS: Eight families with DMJD were studied by whole-exome or targeted sequencing, with detailed clinical and radiological characterization. Patient-derived induced pluripotent stem cells were derived into neural precursor and endothelial cells to study gene expression. RESULTS: All patients showed biallelic mutations in the nonclustered protocadherin-12 (PCDH12) gene. The characteristic clinical presentation included progressive microcephaly, craniofacial dysmorphism, psychomotor disability, epilepsy, and axial hypotonia with variable appendicular spasticity. Brain imaging showed brainstem malformations and with frequent thinned corpus callosum with punctate brain calcifications, reflecting expression of PCDH12 in neural and endothelial cells. These cells showed lack of PCDH12 expression and impaired neurite outgrowth. INTERPRETATION: DMJD patients have biallelic mutations in PCDH12 and lack of protein expression. These patients present with characteristic microcephaly and abnormalities of white matter tracts. Such pathogenic variants predict a poor outcome as a result of brainstem malformation and evidence of white matter tract defects, and should be added to the phenotypic spectrum associated with PCDH12-related conditions. Ann Neurol 2018;84:646-655.

MolProbity: More and better reference data for improved all-atom structure validation.

This paper describes the current update on macromolecular model validation services that are provided at the MolProbity website, emphasizing changes and additions since the previous review in 2010. There have been many infrastructure improvements, including rewrite of previous Java utilities to now use existing or newly written Python utilities in the open-source CCTBX portion of the Phenix software system. This improves long-term maintainability and enhances the thorough integration of MolProbity-style validation within Phenix. There is now a complete MolProbity mirror site at http://molprobity.manchester.ac.uk. GitHub serves our open-source code, reference datasets, and the resulting multi-dimensional distributions that define most validation criteria. Coordinate output after Asn/Gln/His "flip" correction is now more idealized, since the post-refinement step has apparently often been skipped in the past. Two distinct sets of heavy-atom-to-hydrogen distances and accompanying van der Waals radii have been researched and improved in accuracy, one for the electron-cloud-center positions suitable for X-ray crystallography and one for nuclear positions. New validations include messages at input about problem-causing format irregularities, updates of Ramachandran and rotamer criteria from the million quality-filtered residues in a new reference dataset, the CaBLAM Cα-CO virtual-angle analysis of backbone and secondary structure for cryoEM or low-resolution X-ray, and flagging of the very rare cis-nonProline and twisted peptides which have recently been greatly overused. Due to wide application of MolProbity validation and corrections by the research community, in Phenix, and at the worldwide Protein Data Bank, newly deposited structures have continued to improve greatly as measured by MolProbity's unique all-atom clashscore.

Identity-by-descent estimation with population- and pedigree-based imputation in admixed family data.

BACKGROUND: In the past few years, imputation approaches have been mainly used in population-based designs of genome-wide association studies, although both family- and population-based imputation methods have been proposed. With the recent surge of family-based designs, family-based imputation has become more important. Imputation methods for both designs are based on identity-by-descent (IBD) information. Apart from imputation, the use of IBD information is also common for several types of genetic analysis, including pedigree-based linkage analysis. METHODS: We compared the performance of several family- and population-based imputation methods in large pedigrees provided by Genetic Analysis Workshop 19 (GAW19). We also evaluated the performance of a new IBD mapping approach that we propose, which combines IBD information from known pedigrees with information from unrelated individuals. RESULTS: Different combinations of the imputation methods have varied imputation accuracies. Moreover, we showed gains from the use of both known pedigrees and unrelated individuals with our IBD mapping approach over the use of known pedigrees only. CONCLUSIONS: Our results represent accuracies of different combinations of imputation methods that may be useful for data sets similar to the GAW19 pedigree data. Our IBD mapping approach, which uses both known pedigree and unrelated individuals, performed better than classical linkage analysis.

Mechanism of Lysine 48 Selectivity during Polyubiquitin Chain Formation by the Ube2R1/2 Ubiquitin-Conjugating Enzyme.

Lysine selectivity is of critical importance during polyubiquitin chain formation because the identity of the lysine controls the biological outcome. Ubiquitins are covalently linked in polyubiquitin chains through one of seven lysine residues on its surface and the C terminus of adjacent protomers. Lys 48-linked polyubiquitin chains signal for protein degradation; however, the structural basis for Lys 48 selectivity remains largely unknown. The ubiquitin-conjugating enzyme Ube2R1/2 has exquisite specificity for Lys 48, and computational docking of Ube2R1/2 and ubiquitin predicts that Lys 48 is guided to the active site through a key electrostatic interaction between Arg 54 on ubiquitin and Asp 143 on Ube2R1/2. The validity of this interaction was confirmed through biochemical experiments. Since structural examples involving Arg 54 in protein-ubiquitin complexes are exceedingly rare, these results provide additional insight into how ubiquitin-protein complexes can be stabilized. We discuss how these findings relate to how other ubiquitin-conjugating enzymes direct the lysine specificity of polyubiquitin chains.

Molprobity's ultimate rotamer-library distributions for model validation.

Here we describe the updated MolProbity rotamer-library distributions derived from an order-of-magnitude larger and more stringently quality-filtered dataset of about 8000 (vs. 500) protein chains, and we explain the resulting changes and improvements to model validation as seen by users. To include only side-chains with satisfactory justification for their given conformation, we added residue-specific filters for electron-density value and model-to-density fit. The combined new protocol retains a million residues of data, while cleaning up false-positive noise in the multi- χ datapoint distributions. It enables unambiguous characterization of conformational clusters nearly 1000-fold less frequent than the most common ones. We describe examples of local interactions that favor these rare conformations, including the role of authentic covalent bond-angle deviations in enabling presumably strained side-chain conformations. Further, along with favored and outlier, an allowed category (0.3-2.0% occurrence in reference data) has been added, analogous to Ramachandran validation categories. The new rotamer distributions are used for current rotamer validation in MolProbity and PHENIX, and for rotamer choice in PHENIX model-building and refinement. The multi-dimensional χ distributions and Top8000 reference dataset are freely available on GitHub. These rotamers are termed "ultimate" because data sampling and quality are now fully adequate for this task, and also because we believe the future of conformational validation should integrate side-chain with backbone criteria. Proteins 2016; 84:1177-1189. © 2016 Wiley Periodicals, Inc.

The tryptophan metabolism enzyme L-kynureninase is a novel inflammatory factor in psoriasis and other inflammatory diseases.

BACKGROUND: Many human diseases arise from or have pathogenic contributions from a dysregulated immune response. One pathway with immunomodulatory ability is the tryptophan metabolism pathway, which promotes immune suppression through the enzyme indoleamine 2,3-dioxygenase (IDO) and subsequent production of kynurenine. However, in patients with chronic inflammatory skin disease, such as psoriasis and atopic dermatitis (AD), another tryptophan metabolism enzyme downstream of IDO, L-kynureninase (KYNU), is heavily upregulated. The role of KYNU has not been explored in patients with these skin diseases or in general human immunology. OBJECTIVE: We sought to explore the expression and potential immunologic function of the tryptophan metabolism enzyme KYNU in inflammatory skin disease and its potential contribution to general human immunology. METHODS: Psoriatic skin biopsy specimens, as well as normal human skin, blood, and primary cells, were used to investigate the immunologic role of KYNU and tryptophan metabolites. RESULTS: Here we show that KYNU(+) cells, predominantly of myeloid origin, infiltrate psoriatic lesional skin. KYNU expression positively correlates with disease severity and inflammation and is reduced on successful treatment of psoriasis or AD. Tryptophan metabolites downstream of KYNU upregulate several cytokines, chemokines, and cell adhesions. By mining data on several human diseases, we found that in patients with cancer, IDO is preferentially upregulated compared with KYNU, whereas in patients with inflammatory diseases, such as AD, KYNU is preferentially upregulated compared with IDO. CONCLUSION: Our results suggest that tryptophan metabolism might dichotomously modulate immune responses, with KYNU as a switch between immunosuppressive versus inflammatory outcomes. Although tryptophan metabolism is increased in many human diseases, how tryptophan metabolism is proceeding might qualitatively affect the immune response in patients with that disease.

Engineering a genetically encoded competitive inhibitor of the KEAP1-NRF2 interaction via structure-based design and phage display.

In its basal state, KEAP1 binds the transcription factor NRF2 (Kd = 5 nM) and promotes its degradation by ubiquitylation. Changes in the redox environment lead to modification of key cysteines within KEAP1, resulting in NRF2 protein accumulation and the transcription of genes important for restoring the cellular redox state. Using phage display and a computational loop grafting protocol, we engineered a monobody (R1) that is a potent competitive inhibitor of the KEAP1-NRF2 interaction. R1 bound to KEAP1 with a Kd of 300 pM and in human cells freed NRF2 from KEAP1 resulting in activation of the NRF2 promoter. Unlike cysteine-reactive small molecules that lack protein specificity, R1 is a genetically encoded, reversible inhibitor designed specifically for KEAP1. R1 should prove useful for studying the role of the KEAP1-NRF2 interaction in several disease states. The structure-based phage display strategy employed here is a general approach for engineering high-affinity binders that compete with naturally occurring interactions.

Fab-based bispecific antibody formats with robust biophysical properties and biological activity.

A myriad of innovative bispecific antibody (BsAb) platforms have been reported. Most require significant protein engineering to be viable from a development and manufacturing perspective. Single-chain variable fragments (scFvs) and diabodies that consist only of antibody variable domains have been used as building blocks for making BsAbs for decades. The drawback with Fv-only moieties is that they lack the native-like interactions with CH1/CL domains that make antibody Fab regions stable and soluble. Here, we utilize a redesigned Fab interface to explore 2 novel Fab-based BsAbs platforms. The redesigned Fab interface designs limit heavy and light chain mixing when 2 Fabs are co-expressed simultaneously, thus allowing the use of 2 different Fabs within a BsAb construct without the requirement of one or more scFvs. We describe the stability and activity of a HER2×HER2 IgG-Fab BsAb, and compare its biophysical and activity properties with those of an IgG-scFv that utilizes the variable domains of the same parental antibodies. We also generated an EGFR × CD3 tandem Fab protein with a similar format to a tandem scFv (otherwise known as a bispecific T cell engager or BiTE). We show that the Fab-based BsAbs have superior biophysical properties compared to the scFv-based BsAbs. Additionally, the Fab-based BsAbs do not simply recapitulate the activity of their scFv counterparts, but are shown to possess unique biological activity.

Estimating and adjusting for ancestry admixture in statistical methods for relatedness inference, heritability estimation, and association testing.

It is well known that genetic association studies are not robust to population stratification. Two widely used approaches for the detection and correction of population structure are principal component analysis and model-based estimation of ancestry. These methods have been shown to give reliable inference on population structure in unrelated samples. We evaluated these two approaches in Mexican American pedigrees provided by the Genetic Analysis Workshop 18. We also estimated identity-by-descent sharing probabilities and kinship coefficients, with adjustment for ancestry admixture, to confirm documented pedigree relationships as well as to identify cryptic relatedness in the sample. We also estimated the heritability of the first simulated replicate of diastolic blood pressure (DBP). Finally, we performed an association analysis with simulated DBP, comparing the performance of an association method that corrects for population structure but does not account for relatedness to a method that adjusts for both population and pedigree structure. Analyses with simulated DBP were performed with knowledge of the underlying trait model.

Identity-by-descent graphs offer a flexible framework for imputation and both linkage and association analyses.

We demonstrate the flexibility of identity-by-descent (IBD) graphs for genotype imputation and testing relationships between genotype and phenotype. We analyzed chromosome 3 and the first replicate of simulated diastolic blood pressure. IBD graphs were obtained from complete pedigrees and full multipoint marker analysis, facilitating subsequent linkage and other analyses. For rare alleles, pedigree-based imputation using these IBD graphs had a higher call rate than did population-based imputation. Combining the two approaches improved call rates for common alleles. We found it advantageous to incorporate known, rather than estimated, pedigree relationships when testing for association. Replacing missing data with imputed alleles improved association signals as well. Analyses were performed with knowledge of the underlying model.